Non-Gaussian low-frequency noise as a source of qubit decoherence 
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We study decoherence in a qubit with the distance between the two levels affected by random flips 
of bistable fluctuators. For the case of a single fluctuator we evaluate explicitly an exact expression 
for the phase-memory decay in the echo experiment with a resonant ac excitation. The echo signal 
as a function of time shows a sequence of plateaus. The position and the height of the plateaus 
can be used to extract the fluctuator switching rate 7 and its coupling strength v. At small times 
the logarithm of the echo signal is oc t"^. The plateaus disappear when the decoherence is induced 
by many fluctuators. In this case the echo signal depends on the distribution of the fluctuators 
parameters. According to our analysis, the results significantly deviate from those obtained in the 
Gaussian model as soon as w > 7. 
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Introduction: Quantum dynamics of two-level sys- 
tems has recently attracted special attention in connec- 
tion with ideas of quantum information processing. The 
central problem regarding operation of qubits and logical 
gates is maintaining thcphase coherence in the presence 
of a noisy environment jl| . At low temperatures the noise 
is dominated by discrete sources, it is caused by random 
charge exchange between localized states and electrodes 
in the Josephson or semiconductor double quantum- 
dot qubits The charge fluctuations are often mod- 
eled by a set of harmonic oscillators with certain fre- 
quency spectrum In these "spin-boson" models the 
qubit decoherence is determined solely by the pair corre- 
lation function of random forces, Sy{f), that implicitly 
assumes the noise to be Gaussian This assumption 
however does not hold in most practical systems where 
Sx{f) 1// and the processes have extremely broad 
distribution of the relaxation times Q . 

To understand the role of the non-Gaussian statistics, 
we follow and model the environment by a set of two- 
state systems (fluctuators) that randomly switch between 
their states. Their nonequilibrium dynamics can then be 
taken into account explicitly as was done in the anal- 
ysis of coherent quantum transport in the presence of 
1//- noise 0. Recent application of a similar approach 
to qubits demonstrated new features in the decoherence 
that are not reproduced in the Gaussian approxima- 
tion |lOl |. Quantum aspects of non-Markovian kinetics 
were addressed in jllj . 

In the present paper we extend the work |lOl | in two 
directions. Firstly, we evaluate explicitly the phase- 
memory decay in the echo experiment. We find a pro- 
nounced non-Gaussian behavior and explain plateaus 
observed in the time dependence of echo signal 0|- 
Secondly, we consider the case where the interaction 



strengths between the qubit and fluctuators are broadly 
distributed. This distribution strongly modifies the time 
dependence and smears away the plateaus. We suggest 
a recipe for extracting the fluctuators' parameters from 
the measured echo signal. 

It is worth noting that a broad distribution of fluctu- 
ators' switching rates and the coupling strengths makes 
the problem similar to the conventional models of the 
spectral diffusion. It was introduced by Klauder and 
Anderson for the problem of spin resonance. Black 
and Halperin ^3 generalized it 0| to phonon echo and 
saturation of sound attenuation by two-level systems in 
glasses 0. These ideas further developed in 0, 0| 
were applied to the single molecular spectroscopy in dis- 
ordered media ^3 • 

Model: We assume that the qubit is a two-level sys- 
tem (TLS) surrounded by fluctuators - systems with two 
locally stable states. Possible candidates for such fluctu- 
ators in solid state devices are charge traps, see Fig. ^ 
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FIG. 1: Schematic distribution of charged traps located near 
the gate surface and producing oppositely charged images. 



2 



or structural dynamic defects. An occupied trap together 
with its charge image produces a dipole electric field fluc- 
tuating in time due to hops between the trap and the 
gate and acting upon the qubit. If the defects are not 
charged, they can behave as elastic dipoles producing 
time-dependent strains and interacting with the qubit 
via deformational potential. 

A qubit coupled to the environment will be modeled by 
the Hamiltonian H = Tiq+ Hman + "HqF +'Hf, where Hq 
and Tip describe the qubit and the fluctuators separately. 
A completely isolated qubit has two states and is charac- 
terized by the energies of these states and their tunneling 
coupling amplitude. Ti.q is the Hamiltonian of the qubit 
pseudospin cr in a static "magnetic field" , B = {Bx,Bz}. 
Here characterizes the splitting of the energies of the 
two states, and Bx describes their tunneling coupling. 
This Hamiltonian can be diagonalized by rotation in the 
pseudospin space with new z-axis parallel to B. The ro- 
tated Hamiltonian, TCq, is then Hq = —{B/2)(Jz- "Hp 
can be diagonalized by a similar rotation in the fluctu- 
ator's pseudospin space and then split into three parts. 

Hp = Hp^ + Hpi^cnv + Hcnv The first part is just a 
Hamiltonian of an isolated two-level tunneling system, 
Hp"^ = '^i{Ei/2)Tl, where the Pauli matrices r^*^ corre- 
spond to i-th fluctuator. The spacing between the two 
levels, is formed by the diagonal splitting. A;, and 
the tunneling overlap integral, A^, as 



E, = JA^+Af = A,/sin( 



(1) 



The flip-flops of the fluctuators are due to the coupling 
with a thermal bath which we model by a system of equi- 
librium bosons. This applies to phonons, as well as to 
electron- hole pairs in conducting part of the device . 

The interaction, Hqp, between the qubit and the fluc- 
tuators is specified as (cf. with Rcf. 
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Vi = u(ri) cost 



(2) 



Here we assumed for simplicity that the coupling strength 
Vi is determined only by 6i defined by Eq. and the 
distance between the qubit and the zth fluctuator. 

The interaction between the fluctuators and the envi- 
ronment manifests itself through time-dependent random 
fields applied to the qubit. Frequencies of these fields be- 
ing much smaller than the temperature T and qubit split- 
ting _B,the fields can be treated classically: fi*^ — > ^i(t). 
Accordingly, Hqp is the Hamiltonian of the qubit pseu- 
dospin in a random, time-dependent magnetic field X{t) 
formed by independent contributions of surrounding fluc- 
tuators: 

H,F = A'(i)a,, (3) 

i 

The random functions ^i{t) characterize the fluctuators' 
state: ^i{t) instantly switches between ±1/2 at random 



times (random Poissonian process). The switching rates, 
7i, can be calculated in the second order of the pertur- 
bation theory for the fluctuator-phonon/electron interac- 
tion [mill, 



7i = 7o(r) sin^ 9, . 



(4) 



7o is thus the maximal fluctuator switching rate at a 
given temperature, T. For simplicity we assume here that 
the fluctuator spends on average an equal time in each 
state. Although justified only for Ei <C T, this assump- 
tion produces correct temperature dependences |l6j . 

The qubit is manipulated by applied ac "magnetic 
field", F{t) II X, with frequency close to B, so that 
Hman = {^/'2)F{t)a'x- The echo-like manipulation allows 
substantial suppression of the decohcrcncc comparing to 
the "free induction" signal decay . In this case a reso- 
nant ac field is first applied as a pulse rotating the qubit 's 
pseudospin by 7r/2 (7r/2-pulse). After a delay r, a tt- 
pulse inverting the qubit's pseudospin is applied, and the 
echo 7r/2-pulse appears after another delay r. After the 
TT pulse the phase evolves in the reverse direction, thus 
only actual switchings during the time 2r contribute to 
the decoherence, while static fiuctuators do not affect the 
echo signal. 

The external pulses are usually short enough for both 
relaxation and spectral diffusion during each of the pulses 
to be neglected. The echo decay is known to be propor- 
tional to the "phase-memory functional" ^3 



f3{t',T)X{t')dt' (5) 



where (3{t' ,t) = signum(T — t'). The average is calcu- 
lated over the realizations of the random processes ^i(t) 
and random initial states of fluctuators . This averaging 
reflects the experimental procedure where the observable 
signal is an accumulated result of numerous repetitions 
of the same sequence of inputs. Equation (jSJ is obtained 
by analysis of the qubit's density matrix under the per- 
turbation T^man- 

Single fluctuator: Let us start analyzing Eq. ^ with 
the case of a qubit interacting with a single fluctuator. In 
the spirit of Ref. Ila we have obtained the exact solution 
for the echo signal [20|, 
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where /.j = -^/l — {v/2'^Y. This result is essentially non- 
Gaussian. Indeed, assuming that the phase, (pr, obeys 
the Gaussian statistics one would get instead of Eq. Q 
-0(3 = e-<'^')/2 with 



(7) 



Comparing the two expressions we notice that the Gaus- 
sian result (TJ is the weak coupling limit, v <^ 27, of the 
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exact solution ©. However, in the strong coupling case, 
when V > 2j, the exact solution strongly differs from 
Eq. Q. In Fig. 121 both functions are plotted for different 
values of the ratio v/2j. One can see that a.t v > 2j the 




time, yr 

FIG. 2: Echo signal for different values of the ratio w/27 
(shown near the curves), Eq. Dashed lines - calculations 
along the Gaussian assumption, Eq. ||7|l. 

Gaussian assumption strongly underestimates the phase- 
memory functional for t > 2t:/v. Similar conclusion for 
the free induction signal has been recently obtained in 
Ref.ES The reason for the failure of the Gaussian ap- 
proximation in the strong coupling case is similar to the 
well known motional narrowing of spectral lines ■ In- 
deed, if w > 27, then each fluctuator splits the qubit's 
levels rather than broadens them. The qubit just expe- 
riences rare hops between these states; the splitting is 
of the order of v and the typical hopping rate is 7. On 
the contrast, at w <IC 7 the splitting between the levels is 
smeared, the typical decay rate of the echo signal being 
~ v'^ /j. This limit is reproduced within the Gaussian 
assumption. At 7r -C 1, the phase-memory functional 
behaves as 

V'~l-7wVV3, (8) 

regardless of the value of the ratio v/j. This result nat- 
urally holds also in the Gaussian approximation. At 
7T 3> 1 we find 

-ln^«(2^"' ">27; 
^ I w < 27 . 

At u > 27 there appear steps in the time dependence of 
the echo signal shown in Fig. |21 It looks like these have 
been experimentally observed in Ref.0 (see Fig. 3 there). 
At V ^ 7, •y/7/r Eq. ((HJ acquires a simple form 

= 6-2-^^ [1 4- {2-f/v) sinvT] . (9) 

According to Eq. Q, the plateau- like features {dip/dr w 
0) occur at vt ~ 2fc7r and their heights ~ ^-'^-^'^ih 
exponentially decay with the number k. These plateaus 



for 7 ^ w can be understood as follows. The probability 
for the fluctuator to flip during the time 211 jv is very 
small, hence it either does not flip at all or flips only 
once. If it flips during the flrst half of the beating period, 
t < tt/v, the phase of the functional ||SJ) at 2t = 2tt/v 
evolves from to tt. If it flips during the second half of 
the period, tt/v < t < 2tt/v, the phase evolves from tt 
to 27r. After averaging the two contributions will cancel 
each other, which implies that at r close to 2tt/v the 
signal ip is almost insensitive to small variations of t. 

Measuring experimentally the position and the height 
of the flrst plateau, one can determine both the fluctuator 
coupling strength v, and its switching rate 7. For exam- 
ple, the echo signal measured for a Josephson junction 
qubit in Ref.|3 shows a plateau-like feature at r = 3.5 ns 
at the height tp = 0.3. 7 « 27 MHz. If the fluctuator is 
a charge trap near the gates producing a dipole electric 
field, see Fig.^ its coupling strength is v = e^{a- r)/r^. 
Estimating the actual gate-qubit distance r k, 0.5 ^m, we 
obtain a reasonable estimate for the tunnelling distance 
between the charge trap and the gate, a 20 A. 

Many fluctuators: What if the qubit is cou- 
pled to many fluctuators, which are not correlated: 
{S,i(t)(,j(t')) (X 5ij7 The phase memory functional is 
then a product of the partial functionals due to indi- 
vidual fluctuators, \1/ = TT, ip^^^ Following the Holts- 
mark approach Ilfil I2(TI | we approximate Vl^ as 5* ~ 
exp[— ^.(1 — V'^'^)]- Our approach provides accu- 
rate description of the decoherence by fluctuators with 
particular locations as long as the number of active fluc- 
tuators is large. 

We will show now that summation over many fluctu- 
ators with a broad distribution of Vi and 7i may signif- 
icantly change the time-dependence of the echo signal. 
Since the tunneling splitting. A, depends exponentially 
on the distance in real space between the positions of 
the two-state fluctuator, it is reasonable to assume that 
the distribution function P{A,A) = rj/A. In terms of 
E and 6 it implies V{E,6) = rj/smO. Here 77 is pro- 
portional to the number of fluctuators with E < T, 
which are not frozen at a given temperature. Hence, 
?7 oc T. Note also that the contribution of one fluctua- 
tor to the noise spectral density Sx(j-^) is proportional to 
7/(0;^ -1-7^)1 hence the noise spectrum in the frequency 
window 7min <C <C 7o is of the l/u type. 

Due to existence of a finite maximal rate, 70, sum- 
mation over fluctuators with different 7^ does not affect 
much the decoherence at small times. Similarly to Eq. (jH)) 
we find that — In^* cx r'^ for jqt <C 1. This asymptotic 
behavior can however be modified due to a distribution of 
the couphng strengths Vi entering Eq. ((SJ. Let us assume 
that u{r) (X which is the case if the fluctuators act as 
electric dipoles. Uniform distribution of the fluctuators 
over a two-dimensional gate surface corresponds to the 
distribution of the coupling constants V{u) cx which 
coincides with the distribution of the coupling constants 
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in glasses, where two-level systems interact via dipole- 
dipole interaction The phase-memory functional 

can be evaluated by integrating the echo signal ip, 0, 
where 7 and v are related to u and 6 by Eqs. l|2Jl and (0}, 
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T /2 



ir/2 



[1~'4j{u,9)] dd 



(10) 



Here u* is the value of u{r) taken at the average dis- 
tance between the fluctuators having energy splitting 
< T. The integration over u can be extended down to 
zero since {1 ~ xjj) oc at u — > 0. The combination 
C = — (70/u*) In 4* turns out to be an universal function 
of the product jot. At small time C sa (7r/4)(7oT)^. 
A stronger time dependence compared to the single- 
fluctuator case, Eq. is due to strongly-coupled fluctu- 
ators with large v. In the general case, when u{r) oc r~^, 
and the fluctuators are uniformly distributed over an area 
of dimension d, one finds C oc r^+^Z'' for jqt <C 1. Com- 
paring the short-time experimental dependence with this 
prediction one can estimate the ratio d/b, i. c., extract 
an information on the spatial distribution of fluctuators 
and mechanism of their interaction with the qubit. One 
should also keep in mind that energy relaxation described 
by the characteristic time Ti provides an additional echo 
decay as oc e"^"^/-^! , see We neglect this contribution 
to the echo decay. 

In reality, the distribution over u should be cut off at 
some Mmax Corresponding to the minimal possible dis- 
tance between the fluctuator and the gate. Our selection 
of the relevant fluctuators (see above) is valid provided 
that Wmax/70 ^ 1- At 7oT 1 one obtains C en t 0. 

The presented results can be used to probe the sources 
of qubit decoherence by measuring time dependence of 
the echo signal. The plateaus on this dependence signals 
that a major contribution comes from a single fluctuator, 
which parameters can be extracted from the position and 
height of the plateaus using Eq. I^. A smooth time 
dependence of In 5* suggests a combined effect of many 
fluctuators. These conclusions are qualitatively correct 
also for the "free-induction" signal. 

Many fluctuators with a broad distribution of relax- 
ation rates produce noise with l//-typc spectrum. It is 
worth mentioning, see Ref. ^3 for details, that generally 
phase memory decay and 1// noise are dominated by 
different groups of fluctuators. Consequently, in general 
case the fluctuator-induced decoherence cannot be ex- 
pressed through the noise spectrum. This is in contrast to 
the statements that often appear in the literature. Non- 
Gaussian effects due to a single or many fluctuators may 
be responsible for the experimentally observed decoher- 
ence [iii. 
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